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ABSTRACT 

From a Chandra survey of nine interacting galaxy systems the evolution of X-ray 
emission during the merger process has been investigated. It is found that the X-ray 
luminosity peaks ~300 Myr before nuclear coalescence, and then dips, even though 
we know that rapid and increasing activity is still taking place at this time. It is likely 
that this drop in X-ray luminosity is a consequence of outflows breaking out of the 
galactic discs of these systems. In this work it is also shown that, for the systems close 
to the point of nuclear coalescence, Lfir becomes massively enhanced compared to 
the X-ray luminosity of these systems. We suggest that this enhancement may indicate 
a 'top heavy' IMF, with an enhanced fraction of massive stars. 

At a time ~1 Gyr after coalescence, the merger- remnants in our sample are X- 
ray faint when compared to typical mature elliptical galaxies. However, we do see 
evidence that these systems will start to resemble typical elliptical galaxies at a greater 
dynamical age, given the properties of the 3 Gyr system within our sample, supporting 
the idea that halo regeneration will take place within low Lx merger-remnants. 

As a part of this survey, detailed Chandra observations for the double nucleus 
merger system Mkn 266 and the merger-remnant Arp 222 are presented for the first 
time. With the Mkn 266 observation, in contrast to previous studies, we now have good 
spectral information of the individual components part-seen with the ROSAT HRI. 
Additionally, the structure of the emission to the north of the system can clearly be 
distinguished and there is also a suggestion of some extension of X-ray emission to the 
south east of the nuclear region, indicating that this galaxy could just be on the verge 
of large-scale galactic winds breaking out. Within Arp 222 an X-ray luminosity of 1.46 
x 10 40 erg s~ 1 has been detected, this is the lowest value of Lx within our sample. The 
diffuse gas of Arp 222 has been modelled with a temperature of 0.6 keV and, from CO 
observations it has been found to host very little molecular gas, indicating that, from 
current observations, Arp 222 does not resemble a mature elliptical. 

Key words: galaxies: evolution - star: formation - ISM: jets and outflows - stars: 
luminosity function, mass function - galaxies: interactions - X-rays: galaxies 



1 INTRODUCTION 

It is widely believed that very few galaxies exist today that 
have not been formed or shaped in som e way by an interac- 
tion with another galaxy. pToomrel (|l977l ) identified 11 galax- 
ies that exhibit characteristics of on-going mergers, and ar- 
ranged them in a chronological order, illustrating how spi- 
ral galaxies can merge to produce ellipticals. It is now the 
view of many astronomers that this process, whilst not the 
only mechanism to create these systems, plays a vital role in 
the production of elliptical galaxies. The 11 systems within 
the 'Toomre' sequence, alongside other examples of on-going 

* E-mail: njb@star.sr.bham.ac.uk 



mergers, have been studied in great detail over a range of 
wavelengths to try and characterise exactly how these sys- 
tems evolve. When studying these merging galaxies, X-ray 
observations are of particular importance, as they are able 
to probe the dusty nucleus of the system, which can be ob- 
scured at other wavelengths, allowing the nature of the point 
source population to be established. Also imaging of the soft 
X-ray emission permits the diffuse hot gas to be mapped 
out. This gives important information about the hot gaseous 
component associated with the strong starburst, and allows 
galactic-winds outflowing from the system to be observed, 
enabling constraints to be placed on the energetics of these 
outflows. 

X-ray observations were initially carried out with the 
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Einstein Observatory, providing limited spatial resolution. 
This instrument was followed by ROSAT, which greatly im- 
proved the sensitivity of the observations, allowing the X- 
ray properties of these merging galaxies to be probed. A 
study of a sample of int eracting systems was carried out by 
iRead fc Ponmanl (|l998l ) (from here on RP98), where the X- 
ray luminosity and properties of the diffuse gas were investi- 
gated with ROSAT, alongside the point source population. 
In this study it was found that the normalised X-ray lumi- 
nosity of these systems, broadly speaking, followed the nor- 
malised I/fir luminosity, peaking at the time of coalescence. 
It was also found that the young merger-remnant systems 
within the sample exhibited low X-ray luminosities, with 
no indication from the later stage systems in their sample 
that these systems would increase in X-ray luminosity. How- 
ever, it was noted that the ~1 Gyr system, NGC 7252, still 
hosts a large amount of molecular gas, in the form of tails 
and loops, and given more time to evolve, the X-ray proper- 
ties could resemble those of a mature elliptical galaxy at a 
greater dynamical time. The main limitation of this study, 
due to the ROSAT observations, was the inability to disen- 
tangle the point sources from the diffuse gas, particularly at 
larger distances. 

With the next generation of X-ray observatories this 
issue was addressed with an increase in spatial resolution. 
This improvement, provided by Chandra, has allowed the in- 
vestigation of both the diffuse gas and the point source pop- 
ulation of galaxies to be carried out in greater detail. The 
ability to disentangle these two components is vital when in- 
vestigating interacting and merging galaxies, as studies have 
shown that both of these components have different evolu- 
tionary timescales (|Read fc PonmarJ l200ll ). It is therefore 
important to be able to study both of these separately, to 
fully understand how they evolve. 

From observing a selection of interacting galaxy sys- 
tems, at different stages of evolution, the processes involved 
in the merging of galaxies pairs can be characterised. In 
the following sections we will describe the sample we have 
selected, and investigate the behaviour of the X-ray emis- 
sion as the galaxies evolve from two spiral galaxies, through 
to relaxed merger-remnants. Section [2] outlines the selection 
criteria we have used and gives a brief description of each 
system. In section [3] we present the results from new Chan- 
dra observations of Mkn 266 and Arp 222. Correlations and 
evolution of X-ray emission across the whole sample is pre- 
sented in section 3J and discussed in section [5] Conclusions 
are given in section [6] 



2 THE SAMPLE 

To gain a better understanding of galaxy evolution it is im- 
portant to compare similar systems that are undergoing the 
same transformation. By compiling a sample of these galax- 
ies, the evolution of the systems' X-ray properties can be 
investigated. In this paper we have a sample of nine inter- 
acting and post-merger systems, carefully selected to ensure 
that they are representative of galaxies undergoing a major 
merger. The nine systems were selected using the following 
criteria; 

(i) All systems have been observed with Chandra. 



(ii) Systems comprise of, or, originate from, two similar 
mass, gas rich, spiral galaxies. 

(iii) Multi-wavelength information is available for each 
system. 

(iv) The absorbing column is low, maximising the sensi- 
tivity to soft X-ray emission. 

(v) A wide chronological sequence is covered; from de- 
tached pairs to merger-remnants. 

Once the nine systems to study had been selected, the 
issue of merger age had to be addressed. This is one of the 
main problems when working with a chronological study 
such as this, and a number of different methods are required 
to solve this problem. Firstly the point of nuclear coalescence 
was assigned to be at time 0. From this, the time taken un- 
til nuclear coalescence for each pre-merger system can then 
be estimated. These timescales were de rived with a combi- 
nation of N-body simulations, such as iMihos fc Hernquistl 
(1996), and dynamical age estimates, where the length and 
faintness of tidal tails, as well as nuclei separation, were used 
l|Toomre fc Toomrdfl972T) . For post-merger systems, assign- 
ing an age estimate was done by making the assumption 
that the last widespread episode of star formation within 
the system took place at the time of nuclear coalescence. 
Stellar population synthesis models were then used to cal- 
culate the se timescales, therefore g iving a good merger age 
estimate l|Bruzual fc Charlotl 1 19931 ) . In the following sub- 
sections a brief description of the nine merger examples, and 
their X-ray properties, are presented in chronological order. 

2.1 Arp 270 

The earliest example of an interacting system in our sample 
is Arp 270 (also NGC 3395/33 96). This comprises two spiral 
galax ies of comparable mass (Hernandez-Toledo fc Pueraril 
1200 if ), separated by 12 kpc at a distance of 28 Mpc (assuming 
Ho = 75 km s _1 Mpc -1 , and accounting for Virgocentric 
in-fall). These galaxies are connected by an optical bridge 
which is thought to have formed during the systems first 
perigalactic passage which occurred approximately 5xl0 8 
years ago. 

A 20 ks Chandra observation of the system was made i n 
2001 and is discussed in detail in Brassington et al.l (2005), 
the contours of adaptively smoothed 0.3—8.0 keV X-ray con- 
tours overlaid on an optical image are shown in Figure [1] 
From this observation 16 point sources are detected, 7 of 
which are classified as Ultral uminous X-ray Sou rces (ULX's), 
with i x > lxl0 39 erg s _1 (ISoria et al.ll2005l ). The diffuse 
gas emits at a global temperature of ~0.5 keV and shows no 
evidence of hot gaseous outflows as are seen in later stage 
systems (RP98). The galaxy pair, although in a very early 
stage of interaction, already show increased levels of Lpir 
compared to quiescent galaxies, indicating that there is en- 
hanced st ar formation taking place . The numerical simu- 
lations of lMihos fc Hernquistl (|l996l ) suggest that a system 
such as Arp 270 will coalesce in ~650 Myrs. 

2.2 The Mice 

The Mice (also Arp 242, NGC 4676A/B) is another early 
stage merger syste m, lying second in the ev olutionary se- 
quence proposed bv lToomre fc Toomrel (|l972f ). At a distance 
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Figure 1. The early merger system Arp 270. Contours of adap- 
tively smoothed 0.3— S.OkeV X-ray data from Chandra ACIS-S 
overlaid on an optical image from the Palomar 5m telescope. 

of 88 Mpc, the two detached spiral galaxies are again con- 
nected by a tidal bridge and, in addition to this feature, 
exhibit large tidal tails, a consequence of the galax ies inter- 
acting as they have passed each other. iReadl (|2003l ) reports 
on a 30.5 ks Chandra observation of the system, the adap- 
tively smoothed 0.2—10.0 keV X-ray contours overlaid on an 
optical image are shown in Figure [2] 

Five ULXs are detected in association with the galax- 
ies, these sources have been found to be coincident with 
regions of ongoing star formation, both within the nuclei 
of the galaxies and also in the tidal tails. The structure of 
the diffuse X-ray gas in both nuclei suggests that this sys- 
tem is at a more advanced stage of evolution than Arp 270. 
This is indicated by the morphology of the gas, a soft, ther- 
mal plasma, which extends out of the minor axis of both 
galaxies, suggesting that these features are starburst driven 
winds. As in Arp 270, the Mice emits an enhanced level of 
Lfir, again indicating that this system has enhanced star 
formation taking place, this emission is particularly high in 
the nuclei of the galaxies. 

2.3 The Antennae 

The Antennae, NGC 4038/4039, or Arp 244, is probably the 
most famous example of a galaxy pair undergoing a major 
merger, and, lying at a distance of 19 Mpc, is also the near- 
est. A deep integrated Chandra observation of 411 ks has 
been made of this system, enabling the nature of both the 
point source population and the diffuse gas to be investi- 
gated (jFabbiano et aLlbOO^fZezas et al.ll2006t ). Fig |3] shows 
the full band (0.3—6.0 keV) X-ray contours overlaid on an 
optical image. 

Due to the depth of the observation, sources were de- 
tected down to a luminosity of 2—5 xl0 37 erg s _1 . This re- 
sulted in a detection of 120 point sources, 12 of these have 
been confirmed as ULXs. The integrated observation com- 
prises 7 separate pointings, enabling the variability of the 
point sources to be investigated. 

Out of the 12 ULXs, 4 emit below 1 xlO 39 erg s" 1 , 
the lower luminosity threshold for a source to be classified 
as a ULX, in at least one observation. Further, one of the 
sources was observed in only one pointing, indicating that 
it is likely to be a transient, providing further evidence that 
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Figure 2. The Mice, contours of adaptively smoothed 
0.2— 10.0 keV X-ray data from Chandra ACIS-S overlaid on an 
optical image from WFPC2 on board the HST. 

ULXs are a heterogeneous class, comprising of contributions 
from X-ray binary systems, transient sources and also, pos- 
sibly, intermediate mass black holes (IMBHs). 

The quality of the data has also lead to detailed map- 
ping of the diffuse gas, investigating both the temperature 
and metallicity variation of the ISM. 21 separate regions of 
the diffuse gas have been analysed, spectral fitting of these 
regions reveal that there is a variation in temperature of 
the diffuse gas from 0.2 to 0.9 keV and metallicities vary 
from regions of sub-solar abundances to areas of emission 
displaying super-solar abundances, notably in both the nu- 
clear regions and two hotsp ots in the north e rn loo p of the 
disc (Rl and R2 reported in lFabbiano et all (120031 )). 

The morphology of the gas reveals that there are large 
scale diffuse features; two large faint X-ray loops extending 
to the South of the system and a low-surface-brightness halo 
in the region surrounding the stellar discs extending out to 
~18kpc from the nucleus of NGC 4039. The two loops have 
temperatures ranging from 0.29 to 0.34 keV and the low 
surface brightness halo has a nominal temperature of 0.23 
keV. This cooler larger scale emission may be the aftermath 
of a superwind, possibly from the first encounter of the two 
systems, which took place ~2— 5 xlO 8 years ago. The star 
formation rates within the two nuclei are 2.1 Mq yr _1 and 
1.7 Mq yr _1 and in the region where the discs ov erlap has 
been found to be 5.0 M yr _1 (|Mihos et al.ll 1993! ) . 

2.4 Mkn 266 

The next system within this sample, Markarian 266 (also 
NGC 5256), is a double nucleus system with a large gaseous 
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Figure 3. The Antennae, contours of adaptively smoothed 
0.3— 6.0keV X-ray data from Chandra ACIS-S overlaid on an op- 
tical image from the UK Schmidt telescope. 



16:20.0 




13:38:20.019.5 19.0 18.5 18.0 17.5 17.0 16.5 16.0 
Right ascension 



Figure 4. Mkn 266, contours of adaptively smoothed 0.3— 8.0keV 
X-ray data from Chandra ACIS-S overlaid on an optical image 
from WFPC2 on board the HST. 



envelope. Here, a summary of the galaxy properties are pre- 
sented, with the full analysis and results of a recent Chandra 
observation detailed in section I3TT1 

This system lies at a distance of 115 Mpc and has 
an X-ray luminosity of 7.32 xl0 41 erg s _1 . The two nu- 
clei, originating from the progenitors, are thought to com- 
prise a LINER with a powerful starburst component, to the 
north of the system, a nd a Seyfert type 2 to the South 
l|Mazzarella et ai1ll988Tl . In addition to these sources, an 
area of enhanced emission has been detected between these 
nuclei, it is thought that this is caused by the collision of 
the two discs. Prior to a 20 ks Chandra observation of Mkn 
266, made in 2001, this feature had not been seen in X-rays 
before. But, due to Chandra's superior spatial resolution, it 
is now possible to distinguish this feature. This can be seen 
in Figure 31 where the adaptively smoothed, full band X-ray 
contours (0.3— 8.0keV) overlaid on an image from WFPC2, 
on board The Hubble Space Telescope (HST) are shown. 

In addition to the central emission from the colliding 
galaxies, a large X-ray gas cloud can be seen to the North of 
the system. From a ROSAT observation it was proposed 
that this feature is a superwind, driven by the effect of 
the starburst 's supernovae and st e llar winds from the centre 
of the system l|Wang et al.lll997h . iKollatschnv fc Kowatschl 
(1998) argue that an X-ray 'jet', arising from a centrally 
powered superwind, is unlikely due to both the non-radial 
geometry of the emission and also the bright X-ray emis- 
sion of the feature. They suggest that the most plausible 
mechanism to explain the 'jet' is from excitation by hot 
post-shock gas, although they do not conclude where the 
energy to power this feature would arise from. From the 
Chandra observation we speculate that the northern emis- 
sion observed could also arise from star formation taking 
place in a tidal arm that has been stripped from the galaxy 
during an earlier interaction. The exact nature of this emis- 
sion is discussed in detail in section [3.1.41 



2.5 NGC 3256 

The most X-ray luminous system in our sample, NGC 3256, 
is a powerful ultraluminous infrared galaxy (ULIRG), lying 
at a distance of 56 Mpc. This merger system, like Mkn 266, 
has one common gaseou s envelope contain ing the nuclei from 
its two parent galaxies. iLira et all (|2002l ) report on a 28 ks 
Chandra observation made in 2000. Contours of adaptively 
smoothed 0.2—8.0 keV X-ray emission overlaid on an optical 
image of the system are shown in Figure [5] 

The total X-ray luminosity of NGC 3256 is 7.87 x 10 41 
erg s _1 with ~70% of the X-ray luminosity arising from the 
diffuse emission. 14 discrete X-ray sources were detected in 
this system, all of which have been classified as ULXs. Due 
to the high source detection threshold of this observation 
(Lx = 1.4x 10 39 erg s _1 ) fainter point sources were not de- 
tected. Both galaxy nuclei are clearly detected in X-rays. 
The Northe rn nucleus is a site of intense star formation, and 
UV spectra |Lfpari et alj2000h show strong absorption lines, 
implying the presence of massive young stars. The Southern 
nucleus is heavily obscured and appears to be less active then 
the Northern nucleus. It has been suggested that it hosts an 
AGN, although there is no clear evidence of this provided 
in the X-ray data. The soft diffuse emission of the system 
can be described by two thermal components with a harder 
tail. The thermal plasma components exhibit temperatures 
of 0.6 keV and 0.9 keV and the hard component is thought 
to arise from a contribution from the lower luminosity X-ray 
point source population. 

A kinematic study of the system ()English et al.l [20031 ) 
suggests that NGC 3256 is currently experiencing the star- 
burst that just precedes the final core collision and, given 
the close proximity of the two nuclei, it is likely that this 
system has undergone more than one perigalactic approach. 
Assuming that the two tidal tails formed during the last 
closest-encounter, the time that has elapsed since then can 
be estimated. This characteristic timescale was calculated 
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Figure 5. NGC 3256, contours of adaptively smoothed 
0.3— 8.0keV X-ray data from Chandra ACIS-S overlaid on an op- 
tical image from WFPC2 on board the HST. 



to be 500 Myr and it is thought that coalescence will take 
place in ~200 Myr. 



2.6 Arp 220 

The double nuclear system, Arp 220 (IC 4553/4, UGC 
09913) is a prototypical ULIRG, lying at a distance of 76 
Mpc. The two nuclei have a separation of less than 0.5 kpc 
and are just at the point of coalescence. This complicated, 
dusty system has the highest value of Lfir in our sam- 
ple, which arises from the large reservoir of hot gas within 
the system. The mechanism by which this gas is heated is 
still currently under debate. Both the presence of a heav- 
ily shrouded AGN and intense nuclear starbursts have been 
suggested, along with a combination of both these energy 
sources contributing to the heating mechanism. A 60 ks 
Chandra observation of the system was made in 2000, and 
from these data both the nuclear and the extended emission 
have been analysed l|Clements et al.ll2002l ; [McDowel l et al.l 
2003). Figure [6] shows the 0.2—1.0 keV adaptively smoothed 
X-ray emission, overlaid on an optical image. 

The system hosts 4 ULXs, including the nuclei, no other 
point sources are detected. Three distinct regimes of large 
diffuse gas structures are observed; the circumnuclear re- 
gion, a plume region and, further out, diffuse lobe regions. 
The circumnuclear region has both compact, hard nuclear 
emission that resides in the 1 kpc at the centre of the system 
and softer, more extended emission. It has been suggested 
that the hard component of this region arises from a signifi- 
cant point source population, emitting at a luminosity lower 
than ~ 5 x 10 39 erg s _1 , the sensitivity threshold for this 
observation. Although, this could also be a AGN, albeit one 
emitting at a low luminosity. 

The plume regions extend to the northwest and south- 
east of the system with a projected tip to tip length of ~10 
kpc, these features are clearly seen in Figure [6] The spec- 
trum of the plumes include both hot (1—5 keV) and cooler 



(0.25 keV) thermal contributions. It is likely that these re- 
gions are associated with a superwind extending from the 
nuclear region as a consequence of the vigorous star forma- 
tion that is taking place at the centre of the system. Beyond 
the plumes, two large, low surface brightness lobe regions 
have been observed extending from 10—15 kpc on either 
side of the nuclear region. These regions have been found to 
be cooler than the plumes (0.2—1.0 keV), although, due to 
the low number of counts, higher temperature gas residing 

in these regions can not be ruled out . 

It was noted bv lHeckman et all l|l996h that the plumes 
and lobes are "misaligned" by 25°— 30°, they suggest that 
this could be due to a change i n orientation of the syst em as 
the encounter has progressed. iMcDowell et al.l (|2003l ) pro- 
pose that this misalignment is actually a consequence of the 
lobes being produced not by the superwinds in the system, 
but are a product of the merger itself. Tentative evidence to 
support t his scenario is prese nted in an INTEGRAL obser- 
vation bv lColina et al.1 l|2004l ). 

2.7 NGC 7252 

NGC 7252, the first example of a merger-remnant galaxy in 
our sample, is a proto- typical merger galaxy at a distance 
of 63 Mpc. The central part of the system, a single relaxed 
body, displays an r 1 ^ 4 o ptical surface bri ghtness profile typi- 
cal of elliptical galaxies (|Schweizerlll982T ). In addition to this 
relaxed nucleus, the galaxy exhibits complex loops and rip- 
ples, notably two long tidal tails, indicative of the merger 
history of this system. From both U B VI images taken with 
the WCFPT instru ment on the HST (iMiller et all 19971) and 
N-body simulations dHibbard fc M ihos 1995), it has been es- 
timated that nuclear coalescence took place ~1 Gyr ago. 

A 28 ks Chandra observatio n, along with an X MM- 
Newton observation is reported in lNolan et al.l (|2004l ). Fig- 
ure [7] shows the 0.3—7.0 keV adaptively smoothed Chan- 
dra X-ray contours overlaid on an optical image. A total 
of nine ULXs are detected within the optical confines of 
this system, a further 5 sources are also detected but at a 
lower significance (< 3<r). The hot, diffuse gas of the sys- 
tem has been found to be fairly symmetrical and has an 
X-ray luminosity of 2.42xl0 40 erg s _1 . This is low when 
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Figure 7. NGC 7252, contours of adaptively smoothed 
0.3— 7.0keV X-ray data from Chandra ACIS-S overlaid on an op- 
tical image from the LCO 2.5m telescope. 
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compared to luminosities from typical elliptical galaxies (Lx 
~10 41-42 erg s _1 ). The low luminosity of this system is pos- 
sibly due to the young age of the galaxy and, over much 
longer timescales, X-ray halo regeneration may increase the 
mass and hence luminosity of the gas within this system to 
levels seen in typical elliptical galaxies. From spectral mod- 
elling of the XMM-Newton data the X-ray emission from the 
nuclear region is found to emit at 0.72 keV and 0.36 keV. 
There is also a harder contribution that can be fitted with 
a power law, and it is expected that this is due to a lower 
luminosity point source population. 

During the nuclear coalescence of NGC 7252 the star 
formation rate at the centre of the galaxy would have been 
massively enhanced. Now, as the reserve of gas becomes de- 
pleted, the star formation rate has fa llen to one third o f the 
value it would have been at its peak <|Mihos et al.lll993l ), al- 
though its Lfir is still enhanced when compared to normal 
quiescent galaxies. 



2.8 Arp 222 

Arp 222, or 
at a slightly 

NGC 7252 (IGeorgakakis et al 



NGC 7727, is a post 
advanced stage of 
A 



|Crabtree fc Smecker-Hanelll994 ) indicates that these sys- 



200' 



merger galaxy 
evolution than 
CO observation 



tems are very similar in morphology and both host young 
cluster populations, however, Arp 222 seems to have a much 
smaller molecular gas content than NGC 725 2. From a ri- 
band observation ( Rothberg fc Josepbl 12004 ) it has been 
found that the optical surface brightness profile of Arp 222 
follows the de Vaucouleurs r 1//4 law, indicating that violent 
relaxation has taken place since the merger. From analysing 
the discrete structure of the galaxy, plumes extending from 
the nucleus were identified, these features indicate that the 
system has still not fully relaxed into a mature elliptical 
galaxy. 

A 19 ks Chandra observation of the post merger- 
remnant was carried out in 2001, and analysis and results 
from this observation are given in section [3~2l Figure[8]shows 



Figure 8. Arp 222 Contours of adaptively smoothed 0.3— 8.0keV 
X-ray data from Chandra ACIS-S overlaid on an optical image 
from the UK Schmidt telescope. 



the 0.3—8.0 keV adaptively smoothed X-ray emission over- 
laid on an optical image. From this observation 15 point 
sources were detected, two of which are classified as ULXs. 
The X-ray luminosity of the diffuse gas in this system is 6.52 
xl0 39 erg s _1 , much lower than the X-ray luminosity emit- 
ted from NGC 7252. This is likely to be due to the smaller 
gas content of Arp 222, which, as can be seen in Figure [H] 
does not extend to the optical confines of the galaxy. From 
spectral modelling the global temperature of the diffuse gas 
in this system has been found to be 0.60 keV. The Lpir 
value of the system is much lower than that of the previous 
systems and is similar to values seen in elliptical galaxies. 



2.9 NGC 1700 

The final system in our sample, NGC 1700, is a protoellipti- 
cal g alaxy with a best age estimate of ~3 Gyr (|Brown et al.l 
200f|. This galaxy possesses a kinematically distinct core 
and boxy isophotes at larger radii, it also contains two 
symmetrical tidal tails, a good indicator that this system 
formed through the merger of two, comparable mass, spi- 
ral galaxies. It has been suggested that it is the presence of 
these tidal features that causes the 'boxiness' of the galaxy 
IjBrown et aI.ll2000P l. A 42 ks Chandra observation was made 
of the system in 2000, and th e results of this are reported 
in lStatler fc McNamaral (|2002T ). Figure [9] shows the 0.3-0.8 
keV adaptively smoothed X-ray emission overlaid on an op- 
tical image of NGC 1700. 

From this observation, 36 poin t sources are detected , 
6 of which are classified as ULXs (|Diehl fc Statlerl 120071 '). 
The diffuse X-ray gas is well modelled by a single tempera- 
ture thermal plasma at a temperature of 0.43 keV and has 
an X-ray luminosity of 1.47xl0 41 erg s -1 , similar to that 
of a mature elliptical galaxies. The change in morphology 
of the diffuse X-ray gas from the central elliptical region 
to the outer boxy region, can clearly be seen in Figure [9] 
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Figure 9. NGC 1700 Contours of adaptively smoothed 
0.3— 0.8keV X-ray data from Chandra ACIS-S overlaid on an op- 
tical image from the UK Schmidt telescope. 



IStatler fc McNamaral (|2002h suggest that the flattening of 
the isophotes is a consequence of an elliptical-spi ral interac- 
tion, n ot a spiral-spiral merger as suggested bv lBrown et al.l 
(2000). They argue that an interaction involving a preexist- 
ing elliptical with a hot ISM, would lead to the channelling 
of the already hot gas into the systems common potential 
well. This gas, at sufficiently low densities, would then set- 
tle into a rotationally flattened cooling disc, as is observed. 
Currently, neither these two scenarios, nor the suggestion 
that the system could h ave formed from a 3-body interac- 
tion (|Statler et al.l lT996). can be ruled out. 



3.1 Mkn 266 

3.1.1 Overall X-ray Structure 

A 0.3—8.0 keV (from here on referred to as 'full band') Chan- 
dra image was created from the cleaned events file and adap- 
tively smoothed using the CIAO task csmooth which uses a 
smoothing kernel to preserve an approximately constant sig- 
nal to noise ratio across the image, which was constrained 
to be between 2.6 and 4. In Figure fTTJI both the optical im- 
age from the HST with the full band X-ray contours overlaid 
(left) and the 'true colour' image of the galaxy system (right) 
are shown. The 'true colour' image was created by combin- 
ing three separate smoothed images in three energy bands; 
0.3-0.8 keV, 0.8-2.0 keV and 2.0-8.0 keV, using the same 
smoothing scale for each image. These energy bands corre- 
spond to red, green and blue respectively. From these images 
it can be seen that the system is comprised of two separate 
regions of X-ray emission. To the South of the images the 
central emission from the interacting galaxies can be seen. 
This comprises of the two nuclei from the progenitor galax- 
ies, a LINER to the North and a Seyfert 2 to the South, 
contained within a common gaseous envelope and a region 
of enhanced emission between the two nuc lei. This region is 
coinci dent with a radio source reported in iMazzarella et al.l 
(1988). It is likely that this enhanced emission is due to the 
interaction of the two galaxy discs. To the North of the im- 
ages diffuse gas is detected and shows some correlation with 
emission seen in the HST image. The nature of this emission 
is discussed in detail in section [3.1.4l 

From the 'true colour' image it can be seen that the 
north-east nucleus emits in all three energy bands, whilst 
the south-west nuclear emission is not as hard. There is also 
some enhanced X-ray emission in the central region between 
the two nuclei. The X-ray emission surrounding these nuclei 
appears to be soft and diffuse with some suggestion of a 
super-bubble to the south-east of the system. 



3 CHANDRA OBSERVATIONS AND DATA 
ANALYSIS OF MKN 266 AND ARP 222 

Observations of Mkn 266 and Arp 222 were carried out with 
the ACIS-S camera on board the Chandra X-ray Observa- 
tory. Mkn 266 was observed on 2nd November 2001 with a 
total observation time of 19.7 ks and Arp 222 was observed 
on 18th December 2001 with a total observation time of 
19.0 ks. The initial data processing to correct for the motion 
of the spacecraft and apply instrument calibration was car- 
ried out with the Standard Data Processing (SDP) at the 
Chandra X-ray Center (CXC). The data products were then 
analysed using the CXC CIAO software suite (v3.2)Q and 
HEASOFT (v5.3.1). The data were reprocessed, screened 
for bad pixels, and time filtered to remove periods of high 
background (when counts deviated by more than 5<r above 
the mean). This resulted in a corrected exposure time of 
18.5ks for Mkn 266 and 18.8 ks for Arp 222. The full data 
analysis and results from these two observations are detailed 
in the following subsections. 



1 http://asc.harvard.edu/ciao 
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3.1.2 Point Source Spatial and Spectral Analysis 

Discrete X-ray sources were detected using the CIAO tool 
wavdetect. This was run on the full band image, over the 2, 4, 
8, 16 pixel wavelet scales (where pixel width is 0.1"), with a 
significance threshold of 2.5 x 10 -5 , which corresponds to one 
spurious source over a 200 x 200 pixel grid, the size of our 
image. Only 4 sources were detected in this range; the two 
nuclei and two regions in the diffuse northern feature. Both 
the regions in the northern emission contained less than 30 
counts and had large detection regions of r>3.5" and so 
were not defined as point sources. Instead, a spectrum was 
extracted from a region file containing all the northern emis- 
sion. These two detections, along with the extraction and 
background extraction regions for the 2 nuclei, are shown in 
Figure [11] 

The two detected nuclear sources were extracted using 
the source region files created by wavdetect. The size of each 
region was selected to ensure that as many source photons 
as possible were detected whilst minimising contamination 
from nearby sources and background. The background files 
were defined as a source free annulus surrounding and con- 
centric with each source region file, to account for the vari- 
ation of diffuse emission, and to minimise effects related to 
the spatial variation of the CCD response. 
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Figure 10. Left, contours of adaptively smoothed 0.3— 8.0keV X-ray data from Chandra ACIS-S overlaid on an optical image from 
WFPC2 on board the HST. The black box indicates the plate boundary of the HST image. Right shows the 'true colour' image of Mkn 
266. Red corresponds to 0.3-0.8 keV, green to 0.8-2.0 keV and blue to 2.0-8.0 keV. 



The source spectra for the two nuclei were created us- 
ing the CIAO tool psextract and fitted in XSPEC (vll.3.1). 
Due to the low number of counts in this observation the 
Cash statistic (iCashl Il979h was used in preference to \ 2 
when modelling the data. Both sources were well fitted with 
an absorbed thermal model plus an absorbed power law, 
the absorption component was fixed at the value out of our 
Galaxy (1.68 xlO 20 atom cm ). The data were restricted to 
0.3—6.0 keV, as energies below this have calibration uncer- 
tainties, and the spectra presented here do not have signifi- 
cant source flux above 6.0 keV. The parameters from these 
best fit models can be seen in Table [1] where columns 2 and 

3 give the right ascension and declination (J2000), column 

4 the count rate, column 5 the source significance, column 
6 the Galactic value of Nn, column 7 kT, column 8 metal- 
licity, column 9 power law photon index (V) and columns 
10 and 11 give the observed and intrinsic (i.e. corrected for 
absorption) luminosities. 

From this table it can be seen that source A, the LINER, 
exhibits a remarkably flat power law of 0.02±£'ji. It is likely 
that part of the contribution to this component arises from 
unresolved point sources but the very flat slope of this fit 
suggests that some of the hard comp onent is either heavily 
absorbed or dom inated by reflection l|Levenson et al.ll2004l ; 
iMatt et al.ll2000h . By including an absorber for the power 
law component we can allow for this extra absorption, how- 
ever, due to the low number of counts, when including this 
component the fit becomes unconstrained. 

It is likely that the cause of this heavy obscuration is 
du e to the strong starb urst in this region. An observation 
bv lSanders et all (|l986T ) found that Mkn 266 is rich in CO, 
indicating that this system contains a massive, warm reser- 
voir of molecular gas to fuel this starburst, and the dynamic 
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Figure 11. Adaptively smoothed 0.3— 8.0keV X-ray data from 
Chandra ACIS-S of Mkn 266, with wavdetect point source denoted 
with an "X" . Extraction and background region files for source A 
and B are also shown. 

conditions that concentrate this material also serve to fur- 
ther obscure the LINER. The MEKAL component of the fit 
to this source arises from this starburst, which contributes 
5.7xl0 40 erg s _1 to the total intrinsic luminosity of source 
A. 

Source B, the southern nucleus, emits at a temperature 
of 0.88 keV and has a power law slope of 1.13. This thermal 
component is likely to arise from the diffuse emission sur- 
rounding the nucleus, and the value of T is consiste nt with 
the i nterpretation of this nucleus being a Seyfert 2 l|Cappil 
2005). In previous X-ray surveys it has been found that up 
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Table 1. Sources detected in the 0.3—6.0 keV band within the Mkn 266 with a summary of point source spectral fits. Errors on the 
spectral fits parameters arc given as la for 1 interesting parameter from XSPEC. An F denotes that the value has been frozen. 



Source RA Dec. Count Rate Sig. Spectral Fit Luminosity (0.3—6.0 keV) 

N H kT Z T (10 41 erg s" 1 ) 

(xlO" 2 count s" 1 ) (<t) (xlO 20 atom cm" 2 ) (keV) (Z©) Observed Intrinsic 



A 13:38:17.84+48:16:41.2 1.90±0.11 41.8 1.68 F 0.80± '°i 0.59±JJ° 0.02^; 2 f 3.41 ±0.17 3.47 ±0.17 

B 13:38:17.36±48:16:32.5 0.90±0.08 18.7 1.68 F 0.88t ;?o °- 30 F L13 ±a64 °- 98 ±0 - 07 L04 ±0 - 07 



to ~75% of Seyfert 2 objects are heavi ly obscured with JVh 
>10 23 atom cm" 2 (jRisaliti et al.lll999T l. With the data we 
have from this observation the two component model de- 
scribes the spectrum well, and, due to the limited number 
of counts, a more complex model with an additional absorp- 
tion compone nt would over-fit the data. In addition to this, 
ICappil l|2005h have found from a recent survey of nearby 
Seyfert galaxies, that these sources possess the entire range 
of JVh from 10 20 atom cm" 2 to 10 24 atom cm" 2 , fairly contin- 
uously. Indicating that, although the southern nucleus could 
be heavily obscured, it is also possible that the two compo- 
nent model described here is good indication of the proper- 
ties of this source. 

3.1.3 Diffuse Emission Spatial and Spectral Analysis 

From the 'true colour' image of Mkn 266 it is clear that the 
galaxy contains significant amounts of diffuse gas, in both 
the northern feature and also surrounding the central galaxy. 
To investigate the nature of this diffuse emission, spectra 
were extracted from these two separate regions using the 
CI AO tool acisspec, and fitted in XSPEC. Once again, due 
to the low number of counts, the Cash statistic was used 
when modelling the data. 

The diffuse gas contained within the galaxy is not well 
described by a single temperature fit, and is better modelled 
with two temperature components. It is likely that the hot- 
ter gas contributing to this model arises from the enhanced 
emission seen between the two nuclei. To investigate this, the 
extraction region was divided into two separate parts; one, 
smaller region, to probe the impact area, where it is likely 
that the two discs have collided, and a second, larger region, 
covering the rest of the galaxy, but which excludes the inner 
impact region and the northern diffuse region. These, along 
with the extraction region for the northern diffuse emission 
and the associated background regions, are shown in Figure 
1121 The two separate regions, diffuse 1 and diffuse 2, are both 
well described with single component MEKAL fits, exhibit- 
ing the same temperatures that were fitted in the two tem- 
perature model. The northern region is also well described 
with a single component MEKAL fit. The parameters from 
the best fit models for these regions, along with their X-ray 
luminosities, are shown in Table [2] 

From this table it can be seen that the diffuse gas within 
the galaxy system, diffuse 1, exhibits a temperature of 0.52 
keV, a fairly typical temperature for gas in interacting galax- 
ies (RP98). The impact region, diffuse 2, has a higher tem- 
perature of 1.07 keV. This is, within errors, the same tem- 
perature attained from the higher temperature contribution 
from the two component MEKAL fit of the combined spec- 
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Figure 12. Adaptively smoothed 0.3— 8.0keV X-ray data from 
Chandra ACIS-S of Mkn 266, with extraction and background 
region for the diffuse gas shown. The background region files are 
selected to be source free annuli, surrounding and concentric with 
each extraction region file. 

trum for diffuse 1 and diffuse 2, indicating that the higher 
temperature from this fit was a consequence of the gas aris- 
ing from this impact area. 

The northern emission appears spectrally distinct from 
the surrounding diffuse gas and has been found to emit at 
a cooler temperature of 0.30 keV. With this Chandra data 
it has also been found that this region has a much lower 
lumin osity than the one derived in lKollatschnv fc Kowatschl 
(1998), who report on a ROSAT HRI observation of the 
galaxy and find the X-ray luminosity of this region to be 
3.1xl0 41 erg s _1 , as opposed to 6.9xl0 40 erg s -1 , as re- 
ported in this work. However, in iKollatschnv fc Kowatschl 
(1998) this X-ray luminosity is derived from assuming that 
the relative number of HRI counts is directly proportional 
to their share of the integrated X-ray flux, with the total 
luminosity of the system being derived from the ROSAT 
PSPC observation of the system and consequently has large 
uncertainties associated with it. 



3.1.4 Nature of the Northern Emission 

The origin of the diffuse emission to the north of the system 
is still the subject of debate. There have been a number of 
suggestions as to how this feature has been fo rmed. From 
HRI and PSPC observations with ROSAT, IWang et all 
l| 19971 ) suggest that it is an outflow, driven by the mechan- 
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Table 2. Summary of the spectral fits of the diffuse gas regions in Mkn 266. Errors on the spectral fits parameters are given as la for i 
interesting parameter from XSPEC. An F denotes that the value has been frozen. 



Diffuse Region N n kT Z Luminosity (0.3-6.0 keV) 

(xlO 20 atom cm" 2 ) (keV) (Z Q ) (10 40 erg s" 1 ) 

Observed Intrinsic 

Diffuse 1 1.68 F 0.52±npj 0.23±o;°3 9.88 ±0.50 10.92 ±0.50 

Diffuse 2 1.68 F L07 lao9 °- 22 -o.ol K1 ° ±0A9 9 - 83 ±a53 

Diffuse 3 1.68 F 0- 30 tao3 a30F 6 - 10 ±a49 6 - 90 ±a55 



ical energy of the supernov ae and stellar winds in the star - 
burst. A subsequent paper ([Kollatschnv fc Kowatscblll998h . 
investigating both the HRI observation and optical B,V,R- 
images of Mkn 26 6, argued that the scenario proposed by 
IWang et al.l (|l997l ) is unlikely, due to both the high lumi- 
nosity of the 'jet' and also its non-radial geometry. They 
instead suggest that the mechanism that gives rise to the 
'jet' is excitation by hot post-shock gas, although they do 
not conclude where the energy to power this feature would 
arise from. 

A tridimension al spectrophotometric study by 
llshigaki et all (|2000h investigated the Ha, [O in] and 
[S n] emission-lines within Mkn 266. Figure [T3] shows the 
Ha contours from this observation, overlaid on the full 
band, adaptively smoothed, Chandra X-ray emission. As 
can be seen, the enhanced regions of Ha in the northern 
region are coincident with the two X-ray regions detected 
with wavdetect (see Figure lll[) . Some caution should be 
exercised when interpreting this plot due to the uncertainty 
in the astrometry, which we conservatively estimate to be 
2" (1.1 kpc). Even so, the correlation between the X-ray 
and Ha is striking, suggesting that this region to the north 
of the system is a site of star formation, possibly a tidal arm 
that has been stripped from the southern progenitor during 
the merger process. This interpretation is strengthened by 
the fact that both the HST (Figure I10[) and Ha emission 
seem to connect with the d ust lanes around the southern 
nucleus (llshigaki et al.ll2000r i. 

As shown in the previous section, the luminosity arising 
from the diffuse gas within the northern feature is 6.9 x 10 40 
erg s . Given that our spectral fit to the data indicates 
that this emission arises from a thermal component, and 
both the HST and Ha emission indicate that there is star 
formation within this region, the origin of these X-rays could 
come from supernovae. Making some assumptions about the 
geometry of the emitting region of this northern feature the 
thermal energy of the gas can be derived. First the volume, 
V , of the diffuse northern emission has been assumed to be 
an ellipsoid, with symmetry about the longer axis. The fit- 
ted emission measure is equal to r\rfi.V and can be used to 
infer the mean particle number density n e , with the filling 
factor, r), assumed to be 1. This factor represents the frac- 
tion of volume filled by the emitting gas. Although we have 
assumed this to be 1 in our calculations, there is evidence 
from hydrodyn amical simulations to sugges t this value could 
be <2 per cent ^Strickland fc Stevensll2000r ). The mean elec- 
tron density is then used to derive the total gas mass M gas , 
which then leads to the thermal energy Eth of the hot gas. 
From these assumptions Eth has been calculated to be be- 
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Figure 13. Adaptively smoothed 0.3— 8.0keV X- ray data from 
Chan dra ACIS-S of Mkn 266. Ho contours from llshigaki et al.l 
l|2000h are overlaid. 

tween 6.72xl0 54 erg (for 77=0.02) and 3.36xl0 5B erg (for 

77=1). 

By deriving the supernova rate, rsN, for this region, 
the energy a rising from supernovae can be calculated. 
iDavies et all l|2000li calculate r S N for the LINER in Mkn 
266, using the relation derived in ICondon fc Yml (|l990h 

LntIWHT 1 ] ~ 1.3 x 10 23 (7j[GHz])" a8 rsN[yr _1 ], (1) 

and the 20 cm radio continuum (M azzarella et ai1ll98Sf ). 
Where Lnt is the calculated power and v is the frequency 
of the observation. From this they calculate rsN to be 0.45 
SNyr" 1 . 

F or the northern diffuse emission, iMazzarella et all 
( 1988) find a power of 8.9x 10 21 WHz _1 , using this value and 
the a bove equation we fi n d tsn =0.05 SNyr -1 for this region. 
From iMattila fc Meiklel (|200lT ) it has been found that it is 
more appropriate to use Lfir to calculate rsN in starburst 
galaxies but, as we only have Lfir for the whole system, we 
would greatly overestimate the supernova rate, and conse- 
quently, we will use the rate of 0.05 SNyr -1 we have derived 
from the 20 cm continuum da ta. 

From lDavies et al.l l|2000h . models of the star formation 
timescales have been derived for the LINER. And, although 
these models give a fairly poorly constrained age (20-500 
Myr), by making the assumption that this is a proxy for 
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the timescales of star formation taking place in the north- 
ern diffuse emission, the number of supernovae formed in 
this time can be estimated. Given that we also make the 
assumption that a massive star takes ~ 1 x 10 7 yr to evolve 
into a supernova, we calculate the number of SNe formed in 
this time to be 5 x 10 5 — 2.5 x 10 7 . If each supernova re- 
leases ~lx 10 erg, the total thermal energy available will 
be between ~ 5 x 10 56 erg to ~ 2.5 x 10 58 erg, thus demon- 
strating that the observed luminosity in the northern diffuse 
emission could arise from the contribution of star formation 
alone. 

The reason that both IWang et al l <jl997h and 
iKollatschnv fc K owatschl (| 19981 ) prefer an outflow scenario 
to that of tidal stripping to explain this northern feature 
is due to the presence of optical emission in a highly ex- 
cited state within this feature. These line ratios, coupled 
with the electron temperature calculated for this region, 
cannot be explained by thermal collisional ionisation, and 
are likely to arise through ei ther photoionisation or shock 
excita tion. IWang et al.l ()l997h and IKollatschnv fc Kowatschl 
(1998) both suggest that these shocks aris e from a 'jet', 
outflo wing from the central galaxy. However, llshigaki et al.l 
(2000) propose that the northern emission is photoionised 
by radiation from the Seyfert nucleus and therefore, excita- 
tion by shocks is not required to explain the high excitation 
emission lines that are observed. 

The present study adds further evidence to this debate. 
With the higher resolution provided by Chandra, the struc- 
ture of the northern emission can now be resolved. From this 
it can be seen that the northern feature is curved, and the 
X-ray emission appears to be 'clumpy', with the two most 
intense areas of X-ray emission coincident with the Ha emis- 
sion (see Figure I13|l . These are not features which indicate 
a superwind. Furthermore, the morphology of this region 
traces that of the optical emission, which is highly sugges- 
tive that these two features are connected. Consequently, we 
propose that the most likely source of this X-ray emission is 
star formation taking place within a tidal arm that has been 
stripped out from one of the progenitors during the merger 
of the two galaxies. 

3.1.5 The South East Extension 

To the south east of the central galaxy region there is some 
suggestion of X-ray extension, which appears to be coinci- 
dent with filaments seen in the HST image, just beginning 
to break out of the galactic disc (Figure [10)) . Given the low 
count statistics in this observation we cannot extract spectra 
for this region alone. But, from the morphology of the X- 
ray emission, coupled with the Ha emission, there is some 
indication that this region is a site of star formation just 
on the point of break-out from the gaseous envelope of the 
central system. This scenario is preferred to the one sug- 
gested in the case of the northern emission as there is ad- 
ditional evidence of an out flow feature from th e starburst 
region around the LINER l|lshigaki et al.ll200fj| ). This out- 
flow could be beginning to sweep the dust out of the galaxy, 
indicating that this system is in the stage just prior to the 
outbreak of large-scale galactic winds throughout the sys- 
tem. Of course, alternatively, this feature could be a small 
scale version of the northern emission, particularly given its 
coincidence with optical emission seen with the HST. How- 




Figure 15. Adaptively smoothed 0.3— 8.0keV X-ray data from 
Chandra ACIS-S of Arp 222, with point sources detected by 
wavdetect indicated. The D25 ellipse is also shown. 

ever, without spectral information to investigate the nature 
of this object, neither scenario can be ruled out. 

3.2 Arp 222 

3.2.1 Overall X-ray Structure 

A full band adaptively smoothed Chandra image of Arp 222 
was produced using the same tools and techniques as de- 
scribed in section [3. 1.1 1 In Figure [14] both the optical image 
from the UK Schmidt telescope, with the smoothed full band 
X-ray contours overlaid (left), and the 'true colour' image 
(right) of Arp 222 is shown. The true colour image was pro- 
duced using the same methods as described in section 13.1.11 
From these images it can be seen that there is some 
diffuse gas at the centre of the galaxy, but this emission 
does not extend out to the optical confines of the system. 
In addition to the diffuse gas, a number of point sources are 
found through out the galaxy. From the 'true colour' image, 
the white appearance of the nucleus indicates that it emits 
in all three energy bands, whilst the other point sources can 
be seen to be softer. 

3.2.2 Spatial and Spectral Analysis 

X-ray point sources were searched for using the CIAO tool 
wavdetect. This was run on the full band image, over the 2, 
4, 8, 16 pixel wavelet scales (where pixel width is 0.5"), with 
a significance threshold of 2.8xl0~ 6 , which corresponds to 
one spurious source over a 600 x 600 pixel grid, the size of 
our image. 24 sources were detected in this full band range, 
these were then limited to those that lie within the D25 
ellipse of the galaxy, reducing the number to 15 detected 
sources. These regions, along with the D25 ellipse are shown 
in Figure 1151 F rom using the Chan dra Deep Field South 
number counts (jGiacconi et al-lfeoOll ) we estimate 2 to 3 of 
these sources to be background objects. 

Spectra for the point sources were extracted using the 
CIAO tool psextract with region and background region files 
selected in the same way as described in section 15. 1.3 1 Of the 
15 detected sources, only one, source 5, emitted sufficient 
counts to be modelled individually, the other 14 sources were 
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Figure 14. Left, contours of adaptively smoothed 0.3— 8.0keV X-ray data from Chandra ACIS-S overlaid on an optical image from the 
UK Schmidt telescope. Right shows the 'true colour' image of Arp 222. Red corresponds to 0.3—0.8 keV, green to 0.8—2.0 keV and blue 
to 2.0-8.0 keV. 



fitted simultaneously. The spectrum of the diffuse gas was 
extracted with the CIAO tool acisspec, from a region file 
centred on source 5, with a radius of 45". The background 
region file was defined to be a source free annulus surround- 
ing and concentric with the region file. Once again, due to 
the low number of counts, the Cash statistic was used in 
preference to \ 2 - Both the combined sources and source 5 
were well described by single component, power law models. 
And, from this combined fit model, the luminosities for each 
individual source were derived. The diffuse gas was well fit- 
ted by a single component MEKAL model. These spectral 
fits are summarised in Table [3] where column 1 gives the 
source identification number, columns 2 and 3 give the right 
ascension and declination (J2000), column 4 the Galactic 
value of A%, column 5 kT, column 6 metallicity, column 7 
r and columns 8 and 9 give the observed and intrinsic lumi- 
nosities. 

From this table it can be seen that the best-fit spectral 
index for the combined sources is F=1.72, a typical value 
for XRBs in, wha t can be described as, a low/hard state 
ijSoria fc Wul2003h . For the central source, source 5, the best 
fit model has shown that the spectral state is much softer, 
with r=2.15. This value is consistent with it also being an 
XRB, but this time in a high/soft state (|Colbert et alj|2004j ) 
This source, along with source 1, has been found to be a 
ULX. Given the quiescent nature of the galaxy, and that the 
time since the last episode of wide spread star formation ha s 
been shown to be ~1.2 Gyr ago (jGeorgakakis et all [2000) . 
this indicates that these sources are likely to arise from low 
mass X-ray binaries (LMXRBs), not HMXRBs as has been 
seen in younger merger systems within this sample. 

The amount of diffuse gas within the system is smaller 
than in mature elliptical galaxies, and, as mentioned pre- 



viously, does not extend out to the optical confines of the 
system. The temperature of the gas, 0.6 keV, has been found 
to be comparable to that of other interacting systems, but 
the diffuse gas exhibits a lower X-ray luminosity (6.52 xlO 39 
erg s~ ) than the other systems within this sample. The low 
gas content of Arp 222 is further seen from CO observations 
of the galaxy l|Crabtree fc Smecker-Hanelll994h . From this 
Chandra observation it can be seen that this system is X- 
ray faint and does not currently resemble a mature elliptical 
galaxy, although may at a greater dynamical age. 



4 THE EVOLUTION OF MERGING 
GALAXIES 

4.1 Multi-wavelength Evolution Sequence 

To gain a greater understanding of the merger process of 
galaxy pairs, in addition to the X-ray luminosity of each 
system, B-band, K-band and FIR luminosities have been 
obtained to study how each of these luminosity diagnos- 
tics evolve along the chronological sequence. In the case 
of systems where the X-ray luminosity has been obtained 
from the literature, the stated band has been converted into 
the 0.3—6.0 band used in this paper by assuming canonical 
models of the two X-ray components. For the diffuse gas, a 
MEKAL model with a gas temperature of 0.5 keV has been 
assumed, and for the point source population a power law 
model with a photon index of 1.5 has been adopted. 

These luminosities are given in Table|4]where; column 1 
gives the system name, column 2 the distance to the object, 
column 3 the merger age, with nuclear coalescence being 
defined as Myr, column 4 the total (0.3—6.0 keV) intrinsic 
X-ray luminosity from the Chandra observation, column 5 
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Table 3. A summary of the best fit parameters for the diffuse gas and point sources detected in the 0.3—6.0 keV band within Arp 222. 
Errors on the spectral fits parameters arc given as la for 1 interesting parameter from XSPEC. An F denotes that the value has been 
frozen. 



Source 



RA 



Dec. 



(xfO 20 atom cm~ 



Spectral Fit 

kT Z 
(keV) (Z ) 



Luminosity (0.3—6.0 keV) 
(10 39 erg s" 1 ) 



Observed 



Intrinsic 



1 


23:40:00.61 


-12 


17:10.2 


2.75 


F 


1.72F 


0.99 


±0.17 


1.07 


±0.17 


2 


23:39:57.01 


-12 


16:33.0 


2.75 


F 


1.72F 


0.31 


±0.16 


0.32 


±0.16 


3 


23:39:55.95 


-12 


16:47.3 


2.75 


F 


1.72F 


0.19 


±0.07 


0.21 


±0.08 


4 


23:39:54.37 


-12 


16:59.1 


2.75 


F 


1.72F 


0.45 


±0.11 


0.48 


±0.12 


5 


23:39:53.67 


-12 


17:31.6 


2.75 


F 




1.74 


±0.21 


1.97 


±0.24 


6 


23:39:52.84 


-12 


17:36.8 


2.75 


F 


1.72F 


0.37 


±0.11 


0.40 


±0.12 


7 


23:39:52.97 


-12 


18:01.5 


2.75 


F 


1.72F 


0.20 


±0.09 


0.21 


±0.09 


8 


23:39:52.67 


-12 


17:21.1 


2.75 


F 


1.72F 


0.43 


±0.11 


0.46 


±0.12 


9 


23:39:51.67 


-12 


18:19.6 


2.75 


F 


1.72F 


0.52 


±0.14 


0.54 


±0.15 


10 


23:39:51.31 


-12 


18:30.0 


2.75 


F 


1.72F 


0.35 


±0.11 


0.37 


±0.11 


11 


23:39:51.15 


-12 


17:48.8 


2.75 


F 


1.72F 


0.38 


±0.11 


0.42 


±0.13 


12 


23:39:50.57 


-12 


16:53.7 


2.75 


F 


1.72F 


0.28 


±0.09 


0.30 


±0.10 


13 


23:39:51.18 


-12 


15:57.6 


2.75 


F 


1.72F 


0.30 


±0.12 


0.31 


±0.13 


14 


23:39:49.36 


-12 


18:24.4 


2.75 


F 


1.72F 


0.38 


±0.11 


0.40 


±0.12 


15 


23:39:45.63 


-12 


17:41.2 


2.75 


F 


1.72F 


0.21 


±0.12 


0.55 


±0.13 



Diffuse 23:39:53.67 -12:17:31.6 



2.75 F 



0.60 



+0.07 
0.06 



0.08 



+0.06 
0.02 



5.48 ± 0.28 6.52 ± 0.33 



Table 4. The fundamental properties of all the systems within this sample. Columns are explained in the text. 



Galaxy 


Distance 


Merger Age 


£x 


%L dl// 


Log LpiR 


Log L B 


Log L K 


System 






(0.3-6.0 keV) 












(Mpc) 


(Myr) 


(xl0 40 erg s" 1 ) 




(erg s" 1 ) 


(erg s" 1 ) 


(erg s" 1 ) 


Arp 270 


28 


-650 


2.95 


28 


43.63 


43.72 


43.37 


The Mice 


88 


-500 


6.01 


31 


44.12 


44.07 


44.04 


The Antennae 


19 


-400 


6.97 


54 


43.95 


43.98 


44.64 


Mkn 266 


115 


-300 


73.18 


41 


44.75 


44.30 


44.38 


NGC 3256 


56 


-200 


87.70 


80 


45.19 


44.42 


45.22 


Arp 220 


76 





23.70 


14 


45.50 


43.97 


44.78 


NGC 7252 


63 


1000 


6.17 


31 


44.00 


44.40 


44.61 


Arp 222 


23 


1200 


1.46 


45 


>42.37 


43.92 


44.84 


NGC 1700 


54 


3000 


17.58 


83 


42.91 


44.32 


45.15 



gives the percentage of luminosity arising from the diffuse 
gas (%Ldiff), column 6 Lpir, column 7 Lb and column 8 Lk- 
The FIR luminosities ar e calculated using the expres- 
sion l|Devereux fc Ealeslll989t) 



^fir = 3.65 x 10 5 [2.58S , 60Mm + S 100llm ]l 



(2) 



with IRAS 60- and 100-pm fluxes taken from the IRAS Point 

■ 1 — 1| — i 

Source Catalogue (Moshir 1990) . The optical (B) luminosi- 
ties were calculated as in iTullvl Q1988) 



logi B (L Q ) = 12.192 - 0.4B T + 21ogD, 



(3) 



where Bt is the blue apparent magnitude and D is the dis- 
tance in^l£c ; _\^ahMsrf_blue appa rent magnitude were taken 
from Ide Vaucouleurs et alj (Il991 ) (the val ue for The Mice 
was taken from NGC 2000.0 (|Dreverlll988l)). The value s of 
Lk are derived using the relation given in ISeigarl (|2005l ) 

logi K = 11.364 - 0AK T + log(l + Z) + 21ogD, (4) 

where Kt is the K-band apparent magnitude, Z is the 



galaxy redshift and D is the distance in Mpc. Apparent K- 
band magnitudes were taken from the 2MASS survey. 

To calculate the value of %Ldiff for each system, the 
contribution to the luminosity of the diffuse gas that arises 
from unresolved low luminosity point sources had to be es- 
timated and removed. This was done in a number of ways; 
in the case of systems for which we have access to the re- 
duced data (The Mice, Mkn 266 and Arp 222), an additional 
power law component with a photon index of 1.5 was in- 
cluded in the spectral fit of the diffuse emission. It was then 
assumed that the luminosity of this component arises from 
unresolved point sources within the system. Additionally, 
for Mkn 266, it was assumed that the MEKAL component 
of the point source fits arises from the diffuse gas surround- 
ing the LINER and Seyfert. In the cases of Arp 270 and 
NGC 7252, estimates of the contribution of the unresolved 
point so urces to the total diffuse lu mino sity has alrea d y bee n 
made in iBrassington et al.l (|2005l ) and iNolan et all (20041) . 
respectively. 
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For the systems we have taken from the literature, 
the 'Universa l Lum inosity Function' (ULF), derived by 
iGrimm et aD (|2003l ). was used to predict the flux from un- 
resolved low luminosity sources. The predicted XRB lumi- 
nosity function is 

N(> L) = 5.4SFR(L3 8 ' 61 - 21(T ' 61 ), (5) 

where L38 = L/10 38 erg s _1 and the factor 210 arises from 
the upper luminosity cut-off of 2.1xl0 40 erg s . SFR was 
estimated from the population of brighter point sources de- 
tected above the completeness limit from the Chandra ob- 
servations. In the case of The Antennae, point sources were 
detected down to a luminosity of 5xl0 37 erg s _1 and there- 
fore no correction was required for this system. 

As mentioned previously, one of the main difficulties in 
compiling an evolutionary sequence such as this, is assign- 
ing an age to each of these systems. The way in which these 
estimates were made is described in section [2] but another 
important point that was not discussed is that when con- 
structing an evolutionary sample, the absolute timescale of 
the merg er process must be consid ered. From N-body sim- 
ulations dMihos fc Hernquist| [l996) it is clear that from the 
first initial strong interaction between equal-mass mergers, 
through to their nuclear coalescence, takes ~700 Myr. What 
is not so well defined is the amount of time that elapses be- 
tween coalescence and the relaxation of the galaxy into a 
system resembling a mature elliptical galaxy. Within the 
sample presented here, this issue has been addressed by se- 
lecting a greater dynamical range than has previously been 
studied. By doing this it is hoped that the transition be- 
tween young merger remnants to relaxed mature ellipticals 
can be observed, and therefore a more complete picture of 
the merger process can be obtained. 

Although these merger systems are being compared to 
establish a single evolutionary sequence, it should be remem- 
bered that nine individual systems, not nine examples of one 
merger at different stages of its evolution, are being looked 
at. Consequently, although the trends that have been iden- 
tified here should be a good indicator of how X-ray emission 
evolves during the merger process, due to the careful se- 
lection criteria outlined in section [21 it is likely that other 
merger pair systems will exhibit different X-ray properties. 
These variations are likely to arise due to the unique inter- 
action parameters associated with each merger system, as 
well as the variation of gas content and mass of individual 
galaxies within each system. 

With the multi-wavelength luminosities that have been 
collected for each system (Table f3}, the evolution of the 
galaxy properties along the merger process has been in- 
vestigated. The activity levels, a proxy for star formation 
normalised by galaxy mass, is indicated by Lfir/Lk- To 
investigate the variation of Lx, scaled by galaxy mass for 
each system, both the B-band and the K-band luminosi- 
ties, Lx/I/b and Lx/Lk have been used. Both of these have 
been used as normalisation values as Lb can become greatly 
enhanced during periods of star formation due to the pres- 
ence of young stars, indicating that Lk is likely to provide a 
more reliable tracer of galaxy mass. Therefore, by plotting 
Lx against both of these values, how great this effect is can 
be observed. 

These ratios are shown in Figure 1161 where, not only 
the luminosity ratios, but also the percentage of luminos- 



ity arising from diffuse gas (%Ldift ), as a function of merger 
age have been plotted. Lfir/Lk (solid line), Lx/Lb (dot- 
dash line) and Lx/Lk (dashed line) have been normalised 
to the typical spiral galaxy, NGC 2403. Whilst %L dis (dot- 
ted line) is plotted on the right-hand y-axis of the plot, as 
an absolute value, where %L diff for NGC 2403 is 12%. The 
horizontal lines to the right of the plot indicate Lfir/Lk, 
Lx/Lr, Lx/Lk and %Ldig fo r NGC 2434, a typical elliptical 
galaxy (TDiehl fc Statlerll200'it) . 

The system NGC 2403 was selected to represent a typ- 
ical spiral galaxy. This system was chosen due to its close 
proximity (D=3.2 Mpc), the fact that it does not have a 
high leve l of star-forming activity and that it does not host 
an AGN (jSchlegel fc Pannutill2003l ). NGC 2434, was selected 
to represent a typical elliptical galaxy. At a distance of 22 
Mpc, this system is close enough to enable the point source 
population to be disentangled from the hot gas, allowin g 
this diffuse emission to be modelled (|Diehl fc Statlerll2007r i. 
this systems also exhibits the quiescent levels of star forma- 
tion normally associated with such galaxies. Table [5] shows 
the properties for both of these systems, with columns as in 
Table H 

Another factor that must be considered when plotting 
these luminosity ratios is the AGN that is hosted by Mkn 
266. This is the only system within the sample that has 
a confirmed AGN, and as such, it is expected to be signifi- 
cantly more luminous than the other systems. Therefore, the 
values of Lx/Lb, Lx/Lk and %LdiH have been calculated 
for Mkn 266, both including and excluding the contribution 
from the Seyfert. As shown in Table [T] this nucleus is not 
a powerful AGN, instead it is close to the Seyfert-LINER 
borderline, and therefore, the total value of Lx only reduces 
by 14% if it is excluded, with the value of %Ldiff rising to 
48%. When the reduced luminosity is used to derive the lu- 
minosity ratios for Mkn 266, the change in Figure [16] is very 
small, and the overall trends exhibited by these ratios are 
preserved. 

4.2 Arp 220: the Normalisation of the Universal 
Luminosity Function 

IGrimm et all (2003) considered a number of SFR indicators 
for each galaxy within their sample. This was done, as there 
is considerable scatter in the SFR estimates obtained from 
different indicators. To calculate an 'adopted' SFR for each 
galaxy, indicators that deviated significantly from the other 
values were disregarded, and the remaining indicators were 
averaged to give a final value used in subsequent calcula- 
tions. 

In the present work, the SFR was estimated from the 
brighter point sources, detected from the Chandra observa- 
tions. To further investigate the normalisation of the ULF, 
SFR indicators, derived fr om Lfir, using the expression 
l|Rosa-Gonzalez et al.ll2002T ) 

SFRfir = 4.5 x 10- 44 L F iR.(erg s" 1 ), (6) 

have also been used to estimate the luminosity arising from 
the unresolved point sources. This ULF correction, using 
the SFR indicator derived from Lfir, has been calculated 
for the most active system within our sample, Arp 220. From 
this correction, a luminosity of 5.05xl0 41 erg s _1 , has been 
derived for the low luminosity (Lx < 5xl0 39 erg s _1 ) point 
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Figure 16. The evolution of X-ray luminosity in merging galaxies. Shown are Lfir/Lk (solid line), Lx/Lb (dot-dash line), Lx/Lk 
(dashed line) and %-Ldiffi plotted as a function of merger age, where age is defined to be the time of nuclear coalescence. All luminosity 
ratios (Lfir/Lk, Lx/Lb and Lx/Lk) are normalised to the spiral galaxy NGC 2403. %Ldiff is plotted on a linear scale, shown on the 
right-hand y-axis of the plot and is an absolute value. %-Ldiff for NGC 2403 is 12%. The horizontal lines to the right of the plot indicate 
Lfir/Lk, Lx/Lb, Lx/Lk and %Z/diff f° r NGC 2434, a typical elliptical galaxy. 



Table 5. Properties of the typical spiral (NGC 2403) and the typical elliptical (NGC 2434) used in Figure [16] Columns as in Table [I] 



Galaxy 


Distance 


L X 




Log Z-fir 


Log L B 


Log L K 


System 




(0.3-6.0) 












Mpc 


(xl0 40 erg s" 1 ) 




(erg s" 1 ) 


(erg s" 1 ) 


(erg s" 1 ) 


NGC 2403 


3 


0.29 


12 


42.21 


43.22 


43.49 


NGC 2434 


22 


4.82 


76 


41.96 


43.54 


44.49 



sources. This value is almost two times greater than the to- 
tal observed luminosity of Arp 220 (see Table [4) , indicating 
that using Lfir as a SFR in the ULF must greatly overesti- 
mate the point source population for this system. To ensure 
that the value calculated in equation [5] was consistent with 
other observation s, additional Lfir values were obtained 
(Log L F iR=45.86 , |David et all l| 19921 ) and Log L F ir=45.54, 
iLiu fc Kennicuttl ljl995p ). Averaging these values gives Log 
Lfir = 45.73 for Arp 220, similar to the value given in Table 

HI 



To investigate this further, we compare the observed 
population of brighter sources with that predicted using 
the ULF. From the Chandra observation, 4 sources with 
L x > 5x 10 39 erg s" 1 are detected in Arp 220. Using the SFR 
derived from Lfir in equationJS] the bright point source pop- 
ulation is estimated to be ten times larger, with 41 discrete 
sources predicted. This demonstrates that, in the case of Arp 
220, when using Lfir as the SFR indicator, the ULF greatly 
overestimates the whole point source population. Lh q val- 
ues for this galaxy were also used to calculate the SFR of 
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[h] L °S ^ K LFIR 

Figure 17. A plot of the SFR indicated by Lfir and the SFR 
derived from equation [5] with the solid line indicating unity. Arp 
222 and NGC 1700 have not been included in this plot as their 
point source populations are likely to be dominated by LMXRBs. 



the system l|Young et al.lll99fj ; (Colina et alj|2004h but, due 
to the very dusty nature of this object, are subject to large 
extinction effects, leading to greatly reduced SFRs with val- 
ues as little as 0.55, which, given the merger status of this 
system, is unlikely to be a true measure of the SFR. 

To investigate the reliability of using the Lfir SFR indi- 
cator to normalise the ULF, the detected point source pop- 
ulations (with L x < 2.1 x 10 40 erg s~ x ) fr om the Chandra 
observations were used to derive SFRs from equation [5] for 
a number of galaxies in our sample; these values are plot- 
ted against the Lfir derived SFR values in Figure [TT] Only 
seven of the nine systems have been included in this figure, 
as the ULF is only a measure of the HMXRB population, 
and, it is likely that the stellar populations from the two 
older merger-remnants are dominated by LMXRBs. From 
Figure [T71 it can be seen that at high levels of the Lfir SFR 
indicator, the SFR derived from equation[5]is much smaller, 
indicating that, in very active star- forming systems, the pre- 
diction of the point source population, when using the Lfir 
SFR value, does not represent the observed X-ray popula- 
tion well. 

Thi s relationship between Lx and SFR was also investi- 
gated in lGilfanov et all ((2004), where active systems in the 
HDF-N were used to probe the ULF at high values of SFR. 
From this work it was found that, for higher levels of star for- 
mation, the observed X-ray luminosities from these galaxies 
do follow the Lx-SFR relation. However, the calculated SFR 
for all of these systems were derived from 1.4 GHz, L„ SFR 
indicators only, not Lfir SFR indicators. This suggests that 
the discrepancy between Lfir SFR, and the observed values 
derived from the Chandra observations, are a consequence 
of the normalisation indicator chosen, and do not arise from 
non-linearity of the ULF at high values of SFR. 

There are a number of reasons why the SFRs inferred 
from the X-ray observations and Lfir values might differ 
so greatly for systems in the most active phases. It could 
simply be that, for dusty, active systems such as these, Lfir 
is a poor tracer of SFR, becoming greatly enhanced due to 



the large amounts of gas residing in the galaxy. This would 
indicate that it is the Lfir SFR indicator, not the SFR 
from the ULF, that is incorrect. However, this explanation 
is unlikely, as studies have shown that, whilst Lho and Ltjv 
are susceptible to extinction caused by the presence of dust 
in star-forming regions, Lfir does not vary with gas content, 
and scales well with both extinction correc ted Ltjv and radio 
continuum emission in dusty starbursts (|f3uat fc Xul Il996l ; 
lBuatlll992h . 

Alternatively, the discrepancy between the predicted 
and observed values could be a consequence of the age of 
the stellar population. If the starburst has only recently 
taken place, it is plausible that, whilst Lfir has had enough 
time to become massively enhanced, the X-ray point source 
population has not yet evolved into HMXRBs. The conse- 
quence of this would be that the X-ray luminosity arising 
from the point source populati on is lower than the SFR, in- 
dicated by Lfir, would predict. IWilson et ail ((20051 ) recently 
reported on an optical study of Arp 220 which has detected 
two populations of young massive star clusters. The age of 
the youngest of these populations has been found to be 5—10 
Myr, indicating that enough time has elapsed since the last 
starburst for HMXRBs to form (~ 1-10 Myr). It is there- 
fore unlikely that a lag between Lfir and Lx is the reason 
for the difference in the SFR values shown in Figure [T71 

Another reason that the SFR derived from equation [5] 
is lower than expected, when compared to the Lfir SFR, 
could simply be that the X-ray observations are not detect- 
ing all the point sources in very active galaxies due to obscu- 
ration, and, those that are being detected have greatly re- 
duced luminosit ies. To investigate this, the optical obscura- 
tion in Arp 220 (|Shiova et al.ll200lh has been converted into 
an absorbing column density, giving Nn= 2.2xl0 22 atom 
cm -2 , a value which is a ctually somewha t smal ler than the 
modelled value given in IClements et all (|2002[ ) of 3xl0 22 
atom cm" 2 . This demonstr ates that the stated luminosi- 
ties in lClements et all l|2002h are reliable values, and all the 
point sources with Lx> 5xl0 39 erg s _1 (the stated lower lu- 
minosity point source detection threshold) should have been 
detected in this system. 

A final explanation for the difference between the pre- 
dicted and observed SFR indicators for systems close to 
the point of nuclear coalescence, could arise from a physical 
change in the birth of stellar systems. These merging galax- 
ies provide extreme examples of violent star formation, and 
it is therefore possible that the most recent starburst has 
resulted in the birth of a stellar population with a different 
Initial Mass Function (IMF) . 

The IMF, the distribution of masses with which stars 
are formed, has long been a contentious issue, with a variety 
of differe nt models used to e xplain the observed stellar pop- 
ulations (lLarsonlll986l . I 1998). It was initially proposed that 
the IMF was a universal function with a power l aw form, 
regar dless of formation time or local environment (ISalpeterl 
1955). But, it has also been argued th at the IMF varies , and 
was more top-heavy at earlier times l|Rieke et al.lll980l ). 

From the most recent observations it is thought that the 
Salpeter IMF, with a lower mass flattening between 0.5 Mq 
and 1.0 Mq, can be considered a rea sonable approxim ation 
for large regions of starburst galaxies (|Elmegreenl l2005l , with 
some suggestion of small variations with environment, with 
denser and more massive star clusters producing more mas- 



© 0000 RAS, MNRAS 000, 000-000 



The Chandra View of Galaxy Mergers 17 



sive stars compared to intermediate mass stars (|Elmegreenl 
|2004 IShadmehrill2004h . This more top-heavy IMF has the 
effect of generating greater numbers of supernova remnants, 
black holes, and HMXRBs, as well as higher amounts of 
Lfir per unit mass of stars formed, therefore leading to an 
increase in both Lx and Lfir. However, the level of enhance- 
ment of I/fir, relative to the increase in HMXRBs formed, 
is dependent on the shape of the IMF. 

In the galaxy systems close to the point of nuclear co- 
alescence presented in this sample, the X-ray binary pop- 
ulation is dominated by HMXRBs, meaning that Lx will 
scale with the number of massive stars produced, whereas 
Lfir will scale with the main sequence luminosity of those 
stars. As the IMF flattens, Lfir will rise more steeply, a 
consequence of the strong mass to luminosity dependence, 
M oc L 3 ' 5 . Therefore, when using Lfir values to derive the 
SFR for these system, this value will be overestimated. 

If this interpretation is correct, it could explain why 
there appears to be a deficit of X-ray binaries in the Chandra 
observations of the systems close to the point of coalescence. 
In actual fact, what is being seen here is an overestimated 
SFR, derived from the massively enhanced Lfir, as can be 
seen in Figure [TT] 



5 DISCUSSION: X-RAY EVOLUTION 

The first thing to note from Figure \W\ is that Arp 270, 
even though an early stage system, is already exhibiting 
enh anced star format i on ac tivity, as measured by Lfir/Lk 
fsee lRead fc Pon man (200lj)), compared to the typical spiral 
system, NGC 2403. This activity increases up to the point 
of nuclear coalescence, after which there is a steady drop in 
Lfir/Lk until the merger-remnant systems, ~1 Gyr after 
coalescence. From this point the activity value levels off to 
that of a typical elliptical system. 

The evolution of the X-ray luminosity is different to that 
of Lfir/Lk- There is initially a rise, as seen with the star 
formation activity, but this peaks ~300 Myr before coales- 
cence takes place, whilst Lfir/Lk is still increasing. From 
this peak it drops until the young merger-remnants, as is 
the case for Lfir/Lk- But, instead of levelling off, the X- 
ray luminosity once again begins to rise, with the total X-ray 
luminosity of the 3 Gyr system beginning to resemble that 
of a mature elliptical galaxy. 

Both Lx/Lb and Lx/Lk broadly exhibit the same vari- 
ations. One notable exception is the value of Lx/Lb for Mkn 
266. This relatively small value is du e the massive enhance - 
ment of L B from this system's AGN (jRisaliti fc Elvisll2004h , 
which Lk is not susceptible to, therefore indicating that Lk 
is a more reliable tracer of galaxy mass than Lb, and will 
exhibit less scatter. 

The evolution of %Laia also shows an increased value 
for Arp 270, when compared to the typical spiral system. 
This value then steadily rises, up to a point ~200 Myr prior 
to nuclear coalescence. This trend line then drops until the 
young merger-remnant systems, once again rising, with the 
3 Gyr system exhibiting a similar value of %Ldiff to that of 
the typical elliptical, NGC 2434. 



5.1 Previous X-ray Studies 

In RP98 it was found that the X-ray luminosity of the sam- 
ple generally followed the same patterns as the star forma- 
tion activity, with the peak of Lx/Lb being coincident with 
nuclear coalescence. Also in the RP98 study, the sample of 
post-merger systems was limited, with a baseline extending 
out to only 1.5 Gyr after nuclear coalescence. Consequently, 
no rise in the X-ray luminosity of post-merger systems was 
observed, because, as shown in the present study, there is 
a strong suggestion that systems require a much greater re- 
laxation time before they begin to exhibit properties seen in 
mature elliptical galaxies. This idea is further strengthene d 
by a study of post-merger ellipticals (jO'Sullivan et alfcOOll '). 
where a long term trend (~10 Gyr) for Lx/Lb to increase 
with dynamical age was found. 

The reason that the peak in Lx has been found at two 
different merger ages in these studies could be due to the 
different observatories that were used. Chandra has greater 
spatial resolution than ROSAT, and is able to disentan- 
gle background galaxies from the target object with much 
greater accuracy. It is therefore possible that these objects 
were included in the previous ROSAT work as diffuse fea- 
tures, leading to a different peak value of Lx- However, both 
samples include the famous merger system, The Antennae, 
and the system at coalescence, Arp 220, and the X-ray lumi- 
nosities from both of these systems are comparable between 
ROSAT and Chandra. Indicating that this discrepancy does 
not arise from a difference between these two observatories. 

Another difference between RP98 and the current work, 
that could account for the differing merger ages of peak Lx, 
is the selection of systems within each study. In our sample 
we include two systems between The Antennae and Arp 220; 
Mkn 266 and NGC 3256, in RP98 there is only one, NGC 
520. This system is highlighted in RP98 as being X-ray faint, 
not exhibiting the X-ray properties that one would expect, 
given its large Lfir/Lk ratio. They suggest that this galaxy 
does not appear to be on the same evolutionary path as the 
rest of their sample and, possibly, will not evolve into an 
elliptical galaxy. 

A recent study of NGC 520 with Chandra (|Readll2005r ) . 
has indicated that this system comprises one gas-rich and 
one gas-poor galaxy, not two gas rich spirals as has been 
selected within this sample. This lack of gas in the second 
galaxy has resulted in a 'half merger' being induced in this 
system, and hence appears underluminous in X-rays. There- 
fore, it is likely that the peak X-ray luminosity was missed in 
the RP98 sample due to the galaxy selection. In the present 
study, use of the selection criteria outlined in section [2] has 
meant that only systems that are gas rich have been in- 
cluded, ensuring that comparable systems in the merger se- 
quence have been selected. Consequently, the systems pre- 
sented in this sample are more likely to be representative of 
the typical merger evolution. 



5.2 The X-ray Luminosity Peak and the Impact 
of Galactic Winds 

From the X-ray luminosities of the systems within this sam- 
ple it has been seen, for the first time, that both Lx/Lk 
and Lx/Lb peak ~300 Myr before nuclear coalescence takes 
place. This result, given that the normalised star formation 
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rate is still increasing at this time, is initially surprising, as 
one would expect the X-ray luminosity to increase with in- 
creasing Lfir/Lk. However, from looking at the morphology 
of the hot diffuse gas, it is clear that the systems emitting 
very high levels of Lpir/Lk are also experiencing outflows 
from starburst-driven winds. We propose that this relative 
reduction in X-ray luminosity, for these very active systems, 
is a consequence of these large scale diffuse outflows. 

Starburst-driven winds are responsible for the trans- 
port of gas and energy out of star-forming galaxies. The 
energetics of these g alacti c winds were investigated by 
IStrickland fc Stevens! (|2000h . where hydrodynamical simu- 
lations were compared to the observations of the galactic 
wind of M82. From this work it was shown that the major- 
ity of the thermal and kinetic energy of galactic winds is 
in a hot volume-filling component of the gas, which is very 
difficult to probe due to its low emissivity, a consequence of 
the low density of this component. 

Within our sample, the system that represents the peak 
X-ray luminosity is Mkn 266. In section |3~T1 the results from 
the Chandra observation of this system are presented and 
tentative evidence suggesting that this system is just about 
to experience galaxy wide diffuse outflows is found. If this 
interpretation is correct, it seems probable that the lower 
values of L X /L B and L X /L K , for NGC 3256 and Arp 220, 
the most active systems in this sample, are a consequence 
of their extensive galactic winds, which have allowed the 
density, and hence the emissivity, of the hot gas to drop. 



5.3 Halo Regeneration in Low Lx Systems 

Another interesting result from this survey is the increase in 
Lx in the older merger-remnants. In previous studies (e.g. 
RP98) post-merger systems have been found to be X-ray 
faint when compared to elliptical galaxies. In our study the 
merger sequence has been extended to include a 3 Gyr sys- 
tem, and, by doing so, it has been shown that these un- 
derluminous systems appear to increase in L x as they age. 
Given that these merger-remnants have been shown to be 
quiescent, this increase in Lx is not due to any starburst 
activity within the system. Coupling this, with the increase 
in %Ldiff, indicates that diffuse X-ray gas is being produced, 
leading to the creation of X-ray haloes, as obs erved in ma- 
ture elliptical galaxies. lO'Sullivan eiTail (|200lft investigated 
the relationship between Lx/Lb and spectroscopic age in 
post-merger ellipticals and found that there was a long term 
trend (~10 Gyr) for L x to increase with time. The mech- 
anism by which the regeneration of hot gas haloes in these 
galaxies is explained, is one in which an outflowing wind to 
hydrostatic halo phase is driven by a declining SNIa rate. 
lO'Sullivan etafl J200"lh argue that a scenario in which gas, 
driven out during the starburst, infalls onto the existing halo 
is not the dominant mechanism in generating X-ray haloes 
as this mechanism would only take ~1— 2 Gyr and would 
therefore not produce the long-term trend they observe. The 
time baseline from the sample studied in the present paper 
is not sufficient to allow us to discriminate between these 
two possibilities. 




Lo § ( L fir/ L k) 

Figure 18. Plot indicating how Lx/Lk and the two components 
that contribute to Lx/Lk (the point source contribution (L arc ) 
and the diffuse gas contribution (L dif f)) scale with Lfir/Lk. 



5.4 The Behaviour of the X-ray Point Source 
Population 

From the Chandra observations in this study, the behaviour 
of the point source population during the merger process 
has been characterised for the first time. In Figure [18] the 
total luminosity, as well as its two components, the luminos- 
ity arising from the point source population (L src ) and the 
diffuse gas contribution (Ldiff) , have been normalised by Lk 
and plotted against star formation activity (Lfir/Lk). 

From this figure it can be seen that the total Lx/Lk, 
whilst showing a general trend of increasing with Lfir/Lk, 
exhibits a large amount of scatter, peaking at the merger sys- 
tem Mkn 266, and then dropping at high values of Lfir/Lk. 
Decomposing Lx/Lk into separate source and diffuse contri- 
butions to Lx, it can be seen that the scatter arises primar- 
ily from the diffuse component, whilst the drop in Lx/Lk 
at high Lfir/Lk, is a consequence of both Ldiff /Lk and 
Lsi-c/Lk falling steeply. Using the Kendall Rank coefficient, 
it was found that the point source component, L sic /Lk, 
shows a correlation of 2.5a with Lfir/Lk. Fitting a sin- 
gle power law to these data gives a logarithmic index of 
0.93±0.28. 

However, from visual inspection of Figure 1181 a sin- 
gle power law trend does not represent the behaviour ade- 
quately. At high values of Lfir/Lk, the very active systems 
in this sample actually exhibit declining values of Lx/Lk 
as Lfir/Lk increases. In section l4~2l we proposed the pos- 
sibility that these systems close to nuclear coalescence have 
a value of Lfir enhanced relative to that in less extreme 
environments, due to a change in the IMF. Under this hy- 
pothesis, the corresponding points in Figure [18] will have 
been shifted strongly to the right, which could account for 
the negative slope seen in all three curves at high Lfir/Lk. 
The idea that the change in behaviour relates to Lfir/Lk, 
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rather than Lx/Lk, is attractive since, given the very dif- 
ferent mechanisms responsible for generating the source and 
diffuse components of Lx (X-ray binaries and supernova- 
heated gas respectively), it is hard to envisage any process 
which could simultaneously change the behaviour of both. 



6 CONCLUSIONS 

From a sample of nine interacting and merging galaxies, the 
evolution of X-ray emission, ranging from detached spiral 
pairs through to merger-remnant systems, has been inves- 
tigated. As part of this survey, results from the analysis of 
two Chandra observations, Mkn 266 and Arp 222, have also 
been presented. Here we summarise the results from these 
Chandra observations and then draw the main conclusions 
from the survey: 

6.1 Mkn 266 and Arp 222 

• The luminous merger system, Mkn 266, has been 
shown to contain two nuclei, with X-ray luminosities of 
3.47xf0 41 erg s" 1 and 1.04xf0 41 crg s" 1 . In addition, an 
area of enhanced X-ray emission has been detected between 
them. This is coincident with a radio source and is likely a 
consequence of the collision between the two discs from the 
progenitors. A region of diffuse emission is detected to the 
north of the system. It is probable that this arises from a 
spiral arm that has been stripped out of the system during 
the merger. To the south east of the nucleus, a region of 
extended emission has been detected, indicating that this 
system could be on the verge of large-scale galactic winds 
breaking out. This system has the highest Lx/Lk ratio of 
any of the galaxies within our sample. 

• 15 discrete X-ray sources have been detected in Arp 222, 
two of which are classified as ULXs. This merger-remnant 
system has been shown to be X-ray faint when compared to 
both other systems within this sample and the typical ellip- 
tical galaxy NGC 2434. The diffuse gas of Arp 222 has been 
modelled with a temperature of 0.6 keV and, even though 
optical and CO observations are consistent with those of el- 
liptical galaxies, the X-ray luminosity of Arp 222 does not 
resemble that of a mature elliptical. 

6.2 The X-ray Evolution of Merging Galaxies 

• The most striking result from this work is the time at 
which Lx/Lb and Lx/Lk peak. It was previously believed 
that this was coincident with nuclear coalescence, but here 
we find the peak ~300 Myr before this coalescence takes 
place. We suggest that subsequent drop in X-ray emission is 
a consequence of large-scale diffuse outflows breaking out of 
the galactic discs, reducing the hot gas density and allowing 
the escape of energy in kinetic form. 

• This study has also demonstrated that, in the systems 
close to the point of coalescence, Lfir is massively enhanced 
when compared to the X-ray binary luminosity of these sys- 
tems. We suggest here that the high level of Lfir result from 
a change in the IMF in these exceptional starbursts. With 
the production of more massive stars compared to interme- 
diate mass stars in these galaxies leading to larger values of 
Lfir per unit mass of stars formed. 



• At a time ~1 Gyr after coalescence, the merger- 
remnants in our sample are X-ray faint when compared to 
typical X-ray luminosities of mature elliptical galaxies. How- 
ever, we see evidence that these systems will start to resem- 
ble typical elliptical galaxies at a greater dynamical age, 
given the properties of the 3 Gyr system within our sam- 
ple. This supports the idea that halo regeneration will take 
place within low Lx merger-remnants. We caution that, with 
only one older, more relaxed, system within our sample, our 
conclusions on this point are necessarily tentative. To fully 
understand how young merger-remnants evolve into typical 
elliptical systems, the period in which this transformation 
takes place needs to be studied in greater detail. 
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